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ABSTRACT 

This paper investigates the effects of using an approximate state transition matrix in orbit estima- 
tion. The approximate state transition matrix results when higher order geopotential terms in the 
equations of motion are ignored in the formation of the variational equations. Two methods of 
orbit estimation were considered: the differential correction procedure (DC) and the extended 
Kalman filter (EKF). The system used for the study was the Research & Development version of 
the Goddard Trajectory Determination System (R&D GTDS). The effects of the approximation 
were analyzed on a number of orbits. These include orbits of various inclinations and semimajor 
axes. Other parameters studied include geopotential models and DC arc length. 
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In troduction 

The state transition matrix plays an important role in orbit 
determination. It relates perturbations in the state at time t to 
perturbations in the state at epoch. Rice (4) suggests that divergence 
in orbit estimation methods might be linked to the use of an approximate 
state transition matrix. The objective of this project is to study the 
effects of approximating variational equations on orbit estimation 
methods. 

We start with the equation of state: 

A-F(XCOjt) (1) 

which represents n-nonl inear simultaneous equations. An initial state 
vector X(to)=Xo is associated with (1). Tiie state transition matrix is 

described by the matrix differential equation 

# 

(j)= (2) 

Where !3(ta)=I and Fx(X,t) is a matrix of partial derivations of F(X,t) 
evaluated along a particular trajectory satisfying equation (1). 

The force model, F(X,t), used for this study includes perturbations 
involving only gravitational harmonics; other perturbations, such as drag, 
low thrust, etc., have been ignored. Specifically, the force function looks 
like 

■' 1=2. K*0 
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withf(X,t) being the point mass gravitational force caused by the central 
bodyj Oct) the perturbation due to the nonsphericity of the 

i«» 'Cfj) 

central body. The transcendental functions g^(X,t) are extremely complex 

for i‘^3. Based on this force model, equation (2) has the form 

• ^ 

Due to the complexity of g^(X,t), the terms gjj^(X.t) become very cumbersome. + 

The question this study addresses can now be stated as: What is the effect 
on orbit determination methods when M is strictly less than N even though 
the resulting matrix Fx(X,t) is still to be evaluated along a trajectory 
satisfying equation (3)? The main objective for setting MCN in (4) is a 
reduced cost in time and space in programming and evaluating these equations. 

Relationship to Orbit Estimation 

Orbit estimation is the process of solving for the values of a set of 
parameters from the observational model which will minimize the difference 
between a computed and an observed trajectory. The Research Version of the 
Goddard Trajectory Determination System (R&D GTDS) uses two methods of orbit 
estimation: A classical weighted least squares estimator (differential correction 

procedure) and a sequential estimator (Kalman filter). 

The observational model is a nonlinear regression function of the 
state and time; 




+ Baker in Astrodynamics: Applications & Advanced Topics devotes Appendix E 

to "Partial Derivatives of Total Acceleration." 


31 



where n denotes random noise. The system is mxl , m being the number of 
observations. When least squares estimation is considered, the value of X 
which minimizes the weighted sum of the squares of the observational 
residuals is sought. The function to be minimized is called the loss 
function. It has the form: 

^ ( 6 ) 

The initial estimate of the state is Xq. Equation (6) will be minimized 
when 3^“^* Sifice will be nonlinear, Q(X) is first linearized 
by expanding G(X,t) in a truncated Taylor's series about Xq. The linearized 
pt'oblem is now solved and the nonlinear problem is solved recursively via a 
Newton-Raphson iterative scheme to give the minimum difference betv/een computed 
and observed trajectories. This briefly describes the differential correction 
process where a "batch" of m observations are processed simultaneously. The 
state transition matrix is utilized in' the linearization of G(X,t). 

The sequential estimator, or filter, handles the problem from a continuous 
process point of view. Rather than handling the data in batches as in 
differential correction, the filter processes new data immediately upon 
collection to yield an improved estimate of the state. 

In this approach, observations from times tg and t|< are used to determine 
an estimate of the state residual from a reference trajectory X(t|^) and a 
covariance matrix P^. An observation from time t|^+i is added to this set. 
Values of the estimated state at t|(+i , X(t|^+]), and the covariance matrix at 
t|(+i , Pk+1 » to be found. The filter used for this study is the Extended 
Kalman Filter (EKF) as programmed in R&D GTDS. The EKF corrects the reference 
trajectory to the most recent state estimate, which reduces the nonlinearities 
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of the original system and is desirable in real-time sol'jtion., In the EKF, 
the covariance matrix is propagated via the state transition matrix. 


Study Results 

This study has attempted to address the question of approximate state . 
transition matrices by initially investigating a parameter, R, formulated by 

I 

Rice (4). Hr. Rice defines a single parameter to monitor the state transition 
matrix. He presents a statist! caVargument to shov/ that the quantity 

» where (3-jj is an element of the transition matrix (3, can be 
interpreted as a measure of "error growth rate." Ricd' gives P(t)=0P(O)0^ 
as a propagation formula for the covariance matrix, where - 

and states that the square root of the trace of P^t) is commonly used as 
a statistical measure of position errors. Hence, 




The signature of R suggested using the GTDS estimators with several parameters 

to be varied. These included arc length, geopotential modeling of state 

and variational equations, inclination and eccentricity. Being observed 

were the signature of R, the convergence/di vergencd of the estimator, and 

the rate of convergence. 

/ / ' 

As a starting point, three" cases discussed in Mr. Rice's paper were 

compared. Case one used a force model based solely on the point mass force 

for both the state and state transition matrices, which will be>denoted (Jo.Jq)'- 

Case two included the "J 2 " harmonic term in both the equations of state and 

the variational equations, denoted while, ca/e three included the J 2 
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harinonic term 1n the equations of state but only the point mass model in the 
variational equations, comparison of these three cases was 

based on a parameter of the state transition matrix and behavior (in terms of 
convergence/divergence) of the two estimators in R&D GTDS discussed above. 

The first orbit considered was a circular (e-0.0), equatorial (i-0.0°) 
orbit with semi-major axis a=6550.524 km; this orbit will be referred to as 
SATORBl. The EPHEKERIS GENERATION (EPHEM) PROGRAM was used in the calculation 
of the quantity R. EPKEM is used to compute an ephemeris from a given set 
of initial conditions and, optionally, will compute the elements of the 
state transition matrix by numerically integrating the variational equations 
(Eq(4)). Using this option, the quantity R can be printed at any desired 


interval. The first results obtained printed the value of R every 5 minutes 
for the above-mentioned orbit with the modeling of cases 1, 2, and 3. Over 
48 hours, little difference was observed between the corresponding values of 
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R In case 1 (Jq.Jq) and case 2 However, the (J2»Jo^ vastly 

different. While in cases 1 and 2, R grew almost linearly with time, in case 
3, R exhibited an approximately periodic behavior. Repeating the EPHEM runs 
with the same set of initial conditions but a different force model made the 
above comparisons even more striking. In this case, a 4x0 geopotential field 
was used in the equations of state. When generating the partial derivatives, 

4x0, 2x0, and 0x0 force models were used. It is worth repeating at this time that 
the matrix of partial derivatives is always to be evaluated along a trajectory 
of the full equations of motion. Time histories of R for cases (04,04) and 
(04,02) were very similar to those of cases (02,02) and (0 q,0q). Values of R for 
(04,00) followed the same oscillating behavior as (02,0 q). (See Graph 1.) 

These results suggest that when using a model for state with the format of 
eq(3), a simplified force model in the variational equations might be acceptable 
provided that the "O2" geopotential term is explicitly included. 

To test the effects of truncation on batch estimation, a 5-day simulated 
observation file was generated on tape via the GTDS DATASIM program. Range 
and range rate observations were made of SATORBl , where an ephemeris of 
SATORBl was created with J2 included in the state force model. The measurement 
standard deviations for range and range rate were 15 meters and 2 cm/sec, 
respecti vely. With a=6550.625 km, e=. 00012, i = .002°,A =i«>=M=0.0 as an initial 
estimate of the state, DC runs were iriade for cases (J2»>^2) (d2»Jo^ 

hours. In the (J2»J2) case, the DC procedure converged* to the correct solution 
in 4 iterations. The (02»Jo) case diverged. 

These runs were repeated over a 6-hour and 12-hour arc. For the (J2»Jo) 
case, 6 hours was the time at which the parameter R reached its maximum and 


The criteria for converge.ice of the DC are based on the iterative reduction 
of the RMS (square root of the naan square of the observation residuals; 

where are the observation residuals, and m is the number of 


observations. When RMSj+i < RMS- 


the solution is considered converging. 
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hence the time at which R started to decrease in value. In other words, 
for the initial 6 hours, R is monotonically increasing which more accurately 
reflects the "error growth rate" expected in the state transition matrix. 

The period of SATORBl v/as about 90 minutes, so that a 6-hour arc covers about 
4 orbits. When the nc program was used with (J2>Jo) correct over the 6-hour 
arc, it converged in 9 iterations; the 12-hour correction diverged. In other 
words, a short term (where 4xP, P being the period of the orbit, might be a 
guideline for short term) correction might be possible with a point-mass 
force nwdel for the variational equations with a Toss of speed in convergence. 

It has been demonstrated that using an approximate state transition 
matrix can be a detriment to the differential correction process. A logical 
question at this point might be "how much, if any. of an approximation to the 
variational equations can be tolerated by the DC?" The behavior of the parameter 
R when (J 4 ,J 4 ), ( 04 , 02 ) ^nd (Oa,Oo) ephenierides are compared hints that a 
truncation is permissible, provided O2 Is explicitly included in the variational 
equations. To test this hypothesis, simulated observations v/ere made with 
SATORBl elements using a 5x5 geopotential field in the equations of motion. 

Five DC runs were made with Oq, O 2 , O 3 , 0^, and o| models for the variational 
equations and the same initial estimate of the state as mentioned above. The 
DC programs converged in 4 iterations for cases ( 05 , 05 ), (J 5 ^ 0 ^) and (0g,0p. 
Convergence was achieved In 5 iterations for the (o|, 02 ) case. (O^jOq) diverged. 
Table 1 lists RMS values for the last two iterations of these cases as an 
indication of how little is lost when a truncation from O 5 to O 2 for the 
variational equations is used. 


Iteration # 

( 05 , 05 ) 

/ ,5 ,4\ 
( 05 , 04 ; 

(o|,o|) 

(O 5 , 02 ) 

3 

1.1273C19 

1 .0687897 

1.7389176 

2.1896022 

4 

.99626079 

.99626079 

.99631702 

.99647488 

5 




.99626083 


TABLE 1 
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The above runs were all made with a 24-hour arc. The (Ol.Og) <^®5e did 

converge in 9 iterations when used for the short (6 hour) term correction. 

By looking at a typical term of the matrix of partial derivatives, it becomes 

^ k k 

9i^ (X,t) . From Baker, 
the following term is the term added to the 2-body partial derivative when 
forming •^ ^or i=2,3 and K=0: 

To begin with, the term J 2 is three orders of magnitude larger than 3. 

Also, is divided by r"^ , rapidly decreasing the relative magnitude of each 
term for i increasing. In other words, the term J 2 will reflect the vast 
majority of the perturbation due to the asphericity of the earth. Hence, it 
is not at all surprising that the terms , can be truncated when forming 
the partial derivatives without jeopardizing convergence in the DC. 

Other orbits were used to test the relationship between inclination and 
behavior of the DC in the (J2 ,Jq) case. SAT0RB2 had initial elements a=6550.524 km, 
e=.0, 1=30°, n= w=M=0. A DC program was run for a 24-hour arc with an initial 
estimate. of the state as a=6550.624 km, e=. 00012, i=30.002<> and n»o) =M=0. 

Again, the (J2.J2) case converged in four iterations and (J 2 ,Jo) diverged. 

Again, using SAT0RB2 elements with 5x5 geopotential field in the equations of 
motion, a ( 05 , 02 ) DC run over 24 hours will converge in six iterations. These 
results support the suggestion that an approximate state transition matrix might 
be acceptable provided the J 2 potential term is included. 

Several more Inclinations were tried: 1=60O, 1=90°, 1=98°, 1=120°. At 

this point, different results were achieved. The initial orbital elements 
used are listed in Table 2. Simulated observations were made for each orbit. 


Clear that the. J2 term dominates the term ^ 
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a 

e 

i 


XO 

M 

SAT0RB3 

6550.524 

i 

0 

60® 

0 

0 

0 

SAT0RB4 

6550.524 

0 

90® 

0 

0 

0 

SAT0RB5 

6550.524 

0 

98® 

0 

llllllllllll 

0 

SAT0RB6 

6550.524 

0 

120® 

0 

B 

0 


TABLE 2 


For the initial estimate of the state in each DC run, the same error 

was added to the orbital elements: 100 meters added to the semi -major axis, 

eccentricity was increased to .00012, .002® added to the inclination and no 

error added toxi ,io and M. Rapid convergence occured for all (J2>J2) DC 

runs. With computed observations based on an ephemeris of SAT0RB3, the 

(J2,Jo) DC run showed a definite trend toward convergence. After 12 iterations, 

the current state was given as a=6550.257, e'^0(-6), i=60. 00003®, (ar«+M)=720c>. 

DC runs based on simulated observations of SAT0RB4, SAT0RB5, and SAT0RB6 were 
also converging in the (J2>Jo) case, though at a slower rate than SAT0RB3. 

The results are summarized in Table 3. 

Val ue of State 



No. of 
Iterati ons 

3 

e 

1 

SAT0RB4 

20 

6550.534 

,4578x10-4 

89.99754 

SAT0RB5 

27 

6550.525 

.2398x10-5 

98.00005 

SAT0RB6 

■ 1 

24 

6550.252 

.6327x10-5 

120.0 


TABLE 3 
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Table 4 compares the RMS value for various iterations for SAT0RB3, 4, 
5, and 6 (J2.Jo)* serves as a monitor for the rates of convergence. 


RMS Values (Og.Jo^ 


Iteration # 

SAT0RB3 

SAT0RB4 

SAT0RB5 

SAT0RB6 

1 

2084.7695 

2261.6206 

1766.1576 

2307.2376 

6 

74.458092 

1175.0753 

742.93831 

103.57066 

12 

57.808147 

521.06515 

260.84572 

112.82074 


TABLE 4 


As one might expect, since cos 5^ =\cos^| and sin^ = sin^ 

the RMS values are most similar for 1=60° and 1=120°, (SAT0RB3 and SAT0RB6). 

In these two cases, the first 12 Iterations alternated between converging and 
diverging, with large decreases of RMS value in convergent Iterations. 

(This accounts for RMS ^2 ^ in SAT0RB6.) 

In order to examine the sensitivity to Inclination, it is helpful to 
look at first order perturbations. In Methods of Orbit Determination . Escobal 
devotes a chapter to "Secular Perturbations," where the term secular describes 
variations "associated with a steady nonscillatory , continuous drift of an 
element from the adopted epoch value."* He represents the perturbing potential 
as Asl-V where 3T Is the potential due to an aspherical earth and V is the 
potential of a spherical earth. He segregates from A those terms which 
will contribute secular variations In the elements and arrives at 

= (7) 

where K^nm^a^. Note that this is a first order expression in J 2 and for the 
sake of this analysis, the , 1 3, terms have been neglected. Little Is lost 

* Escobal, P.R., 1965, p. 362. 
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by neglecting J-j , 1 ^ 3 as the J 3 term is approximately 10 orders of magnitude 
smaller than the J 2 term and the relative magnitude of Ji and J 2 becomes even 
more drastic for 3. Expression (7) is then averaged over one revolution, 
resulting in: 



Using this as the perturbing function due to J 2 , it is easy to see that the 
secular effect of J 2 is eliminated in the equations of motion v/hen i=54.70 
(since 1/3-1/2 sin^ 54.7O=0). In other words, at this specific inclination, 
the satellite, in a secular sense, perceives the earth as approximately (i.e., 
to first order J 2 ) spherical. 

With the aid of the above model for the perturbing function, Escobal 
develops the following equations representing the gradual drift of the classical 
elements from their adopted epoch values. Note that only-fx,'tO and M experience 
this drift and a, e and i are taken to be constant. (It might be worthwhile to 
state again that this is only a first order secular perturbation theory.) 

Anomal i sti c mean moti on : 

ri = YVt Ivi 

Mean Anomaly: 

M - (•t-to') 

Longitude of the ascending Node: 

_fL s= ^6 - n (t'-fco) 

Argument of Perigee: 

ox)* Uc'*' 
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Between these critical Inclinations of 54. 7o and 125.3°, the rate of the 
secular variations of the elements is smaller than outside of this region. 
This accounts for convergence in the DC procedure with niodeling 

for orbits with inclinations betv/een 54° and 125°. 

To conclude the differential correction section of this study, two more 
orbits were considered. These orbits were elliptical with a greater 
semi -major axis. Simulated observations were made of these orbits. The 
initial states are given in Table 5. 



a 

e 

i 


X*.) 

M 

S.AT0RB7 

7278.360 

.1 

0 

0 

B 

0 

SAT0RB8 

9357.89143 

.3 

0 

D 

B 

0 


TABLE 5 


A larger value for a will decrease the effect of J2 which is readily seen in 
equation (8). However, the effect of J 2 1 ^°^ absent from these two orbits 
and the graph of the parameter R with (J2 ,Jq) modeling suggests that the DC 
procedure will have trouble converging over a 24 hour arc, which it indeed does. 
But the period of SAT0RB7 and 8 is increased to 103 minutes and 150 minutes, 
respectively. Because of its period, it is not surprising that SAT0RB8 converges 
in 11 iterations over a 12-hour span with the (J2 ,Jq) modeling. With SAT0RB7, 
P=103 minutes so that 7 hours (<v4xP) should be a reasonable time arc in the 
(02,Jo) When a 6-hour arc is used, the (02>'^o) DC converges to SAT0RB7 
elements in 11 iterations with rapidly decreasing RMS values. With a 12-hour 
arc, convergence was still not achieved after 30 iterations. 

The results obtained using the FILTER as an orbit estimator are more 
difficult to examine than the results from the DC. As input to the FILTER 
program, the user supplies an initial estimate of the state along with an 
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Initial estimate of the covariance matrix.* The a priori covariance 
matrix contains the state standard deviation and correlations, hence 
points the filter in the right direction. 

For this study, the observational residuals were used to monitor 
filter performance, with decreasing residuals within a pass and modestly 
larger residuals appearing after a data gap indicating convergence. The 
arc length used in the filter portion of this study was 18 hours. 

When testing the effect of approximating the state transition matrix 
in the filter, only the SATORBl orbit was used. With i=0.0° and e=0.0, 
this orbit is particularly susceptible to perturbations due to the earth's 
J 2 nonsphericity. As will be demonstrated below, the filter has an added 
dimension of sensitivity, that being the a priori covariance matrix. Because 
of this and time constraints the use of the filter was restricted to this 
orbit. 

The first offset imposed on the state was the same as that used for 
the DC:4 a=100m,Ae=. 00012, ai = . 002° , V =ii=M=0.0. With this offset, 
the cartesian elements at to are x=6549.8379 km, y=z=0.0 km, x=0 km/sec, 
y=7. 801531422 km/sec and 2 =-. 0002723 km/sec. At tg. the true cartesian elements 
are x=6550.524 km, y=z=0.0 km, x=z=0.0 km/sec and y=7. 8006548 km/sec. Four 
different a priori covariance matrices were tried in the filter with this 
initial state estimate. All four covariance matrices were diagonal , implying 
there was no correlation among the errors in the state estimate. The first 
covariance matrix exactly reflected the errors in the state: 

* The DC procedure has an a priori covariance matrix default value of 
infinite magnitude so that its inverse is the null matrix. In the DC 
procedure it is this inverse which is an additive term to the loss function 
(eq(6)); however, it has been omitted in eq(6) as it is the null matrix. 
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(T* x=.68 km,<ry=o -'2 = .1 X x 10-^ krn/sec, cr^=.87 x 10"3 km/sec 

and 0~j=.27 x 10"^ kni/sec. 

fhe results obtained for (J2>J2) ^rid (02>Jo) were very similar, both cases 
converging. Table 6 below lists, for the last 5 passes, the largest 
residual in meters within a pass for both cases: 

Pass Number 



19 

20 

21 

22 

23 

24 

(J2.J0) fnax 10~c! 

23 

20 

29 

21 

19 

31 

( 02 , 02 ) max |0“Cl 

31 

43 

22 

1 

21 

30 


TABLE 6 

These results change greatly v/hen a different a priori covariance 
matrix is used. With km,(J~2=.01 km and ^=<5^=0"^=. 1 x lO"'^ km/sec, 

the (J2 >Jq) f-'ise failed to converge. This covariance matrix fails to 
recognize the error in the y component, Ay=. 00087 km/sec, so that the actual error 
in y is 87 times larger than is reflected in the standard deviation associated 
with it: IT y=.l x 10'^ km/sec. The largest residual in the last 6 passes is 
listed below. Note that the (J2,J2) ^^res much better than (O 2 .O 0 ) and 
is considered converging. 

Pass Number 



19 

20 

21 

22 

23 

24 

(02,0q) max ] 0-Cl 

997 

1151 

1104 

1094 

1250 

982 

( 02 , 02 ) max )0-Cl 

35 

72 

21 

58 

45 

81 


TABLE 7 

If the standard deviation of y is increased to 0” y=. 00004 so that y is 
approximately 20 times larger than is reflected in the standard deviation. 
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the (J2»'^0) responds a little better. However, it is not considered 
converging. The residuals during the first 9 hours suggest the filter 
has a handle on the correct state. The next 9 hours shows increasing 
residuals as the filter drifts away from the possible steady state. 

Lastly, the standard deviation of y was increased again to 0"y=. 00015 or 5 times 
smaller than the error on y. Here the (J 2 »do) esse did as well as the 
(J2,J2) case. Table 8 summarizes the results for <3"y=.0004 and 
<Ty=.00015, with rj^=(Ty=.l,<5‘z=.01,!rx=(rz=.l x lo“^. 


Pass Number 



19 

20 

21 

22 

23 

24 

=.00004 km/sec 
(J?. jJq) 

190 

226 

2.30 

220 

242 

227 

(J2,J2) 1 0-Cl 

28 

39 

22 

28 

12 

41 

=.00015 km/sec 

(02 »>lo) 1 

32 

23 

31 

33 

29 

47 


28 

37 

23 

28 

14 

38 


TABLE 8 


Although these results are far from conclusive, some inferences can 
be drawn from them. It has been demonstrated that for this case the filter 
responds well with an approximate state transition matrix provided the a 
priori covariance matrix reflects the state errors within five standard deviations. 
The test cases used for this study suggest that when the error is between 20 
and 90 tioes larger than the standard deviation, the (J2,0o) fails, yet 
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the (J2»'j2) '*'0 >"each a steady state solution. This suggests 

that when the a priori covariance matrix is considered an accurate indication 
of the state error, truncated variational equations might not harm the per- 
formance of the filter. On the other hand, v/hen the initial covariance matrix 
somewhat inaccurately reflects the state error, the full partial derivatives 
are needed to help steer the filter toward steady state. (The term "somewhat" 
is used here as a precaution; a totally inaccurate a priori covariance matrix 
can easily cause a (J 2 .J 2 ) filter case to diverge.) This sensitivity to the 
initial covariance matrix makes it difficult to draw conclusions regarding 
the filter's performance as a function of the state transition matrix. 

Conclusions 

This paper has attempted to evaluate the effects of an approximate 
state transition matrix on the differential correction procedure and the 
filter procedure as used for orbit estimation. The DC results fall into 
four categories: the effects due to (1) extent of the approximation, 

(2) orbital inclination, (3) length of time arc, and (4) orbital 
eccentricity. Coinciding with these categories is the behavior of the 
parameter • R can be used to "predict" convergence/divergence 

in the DC and its behavior suggested these four categories as meaningful 
avenues to investigate. When using a force model for the state which 
includes the j|JJ harmonic term, it has been shown that it is a safe practice 
to approximate the variational equations provided that the J 2 term is 
included. When this approximation is made, the DC process will still 
converge as it would with the full variational equations with only a 
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negligible loss of speed. A total truncation of the harmonic terms in 
the variational equations will, in general, cause divergence in the DC. 

One exception to this is orbits with inclinations between 54.7° and 125.3° 

In this range of inclination, the effect of the nonsphericity of the earth 
is minimized. Here (J2 >Jq) cases will converge but so slowly that 
the truncation might, in practice, be undesirable 

The oscillatory signature of R when the most drastic truncation is 
made suggests the possibility of short term differential corrections. The 
maximum value of R tends to occur after four periods of the orbit. Hence an 
arc length of four times the orbital period becomes a reasonable guideline 
for short term corrections. In this arc length, convergence is achieved but 
speed of convergence again becomes the trade-off for the truncation. As the 
orbital period lengthens with greater semi-major axis, so does the time 
arc over which the (J 2 .J 0 ) converges in the DC procedure. 

With the filter as an estimator, less conclusive results are found. 

The effect of a truncation in the variational equations on filter performance 
was highly correlated to the initial covariance matrix. When the. a priori 
covariance matrix was a good indication (within 5 standard deviations) of 
the actual error imposed on the state, little difference was seen in the convergent 
behavior of the filter for the (J2»J2) (d2>'^o) cases. However, when the 

accuracy of the a priori covariance matrix is relaxed, the (J2»Jo) showed 
divergence in the cases tested. Furthermore, when the accuracy of the a priori 
covariance matrix is completely lost, neither the (J2»Jo) (J2»'^2^ 

will converge. 
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The results obtained in the filter section of this paper leave another 
set of questions open. It was assumed that using the low altitude, circular 
equatorial orbit where the perturbations due to the nonsphericity of the earth 
are most pronounced would be a good orbit to test the effects of truncated 
variational equations in the filter. This appears to be a valid assumption 
in view of the analysis mentioned above based on equations (7) and (8). 
However, other sets of orbital elements could be tested in order to help 
determine the limits of accuracy needed in the a priori covariance matrix 
when using an approximate state transition matrix. Also, this study was 
restricted to variations inO”i^ . Certainly many other variations could be 
tested, although this starts to drift away from the original intent of the 
study. Also, is there a level of state noise which might be used to help 
compensate for the use of an approximate state transition matrix? 

In general, this study could be expanded to nonpotential accelerations 
such as drag and solar radiation pressure. What, then, would be the effect 
of including a nonpotential acceleration in the equations of motion but 
excluding it in the variational equations? 

Lastly, the stability properties of the state transition matrix is 
a question of interest. Does the solution to the variational equations 
exhibit one type of stability for the (J2,J2) modeling which is different 
from (J2,Jo) modeling? 

Although many questions remain open, it is hoped that this study 
sheds some light on the appropriateness of state transition matrix 
approximations. 
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